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SIMMARY 


The minimum critical Reynolds numbers for the similar solutions of 
the conqpressible laminar boundary layer computed by Cohen and Reshotko 
and also for the Falkner and Skan solutions as recomputed by Smith have 
been calculated by Lin's rapid approximate method for two-dimensional 
disturbances. These results enable the stability of the compressible 
laminar boundary layer with heat transfer and pressure gradient to be 
easily estimated after the behavior of the boundary layer Viar been com- 
puted by the approximate method of Cohen and Reshotko. 

The previously reported unusual result (NACA Technical Note 4057) 
that a highly cooled stagnation point flow is more unstable t.)mn a highly 
cooled flat-plate flow is again encountered. Moreover, this result is 
found to be part of the more general result that a favorable pressure 
gradient is destabilizing for very cool walls when the Mach number is 
less than that for complete stability. The minimum critical Reynolds 
numbers for these wall temperature ratios cure, however, all larger than 
any value of the laminar-boundary-layer Reynolds number likely to be 
encountered. For Mach numbers greater than those for which complete 
stability occurs a favorable pressure gradient is stabilizing, even for 
very cool walls. 


INTRODUCTION 


In reference 1 a useful method for calculating the compressible 
laminar boundary layer with heat transfer and arbitrary pressure gradient 
is presented. This method is based on the similar solutions of the lami- 
nar boundary-layer equations obtained in reference 2. 

Because of the importance of the problem of transition from i>m^Ar 
to turbulent flow, it is often desirable to have an estimate of the sta- 
bility of the laminar boundary layer. In order to obtain such an estimate 
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easily, the minimum critical Reynolds numbers for the similar solutions 
presented in references 2 and 3 have "been calculated for the Mach number 
range between 0 and 2.8 by the rapid approximate method of reference 4. 
The results are presented in tables and charts so that, after a calcula- 
tion of the laminar boundary has been made by the method of reference 1, 
the distribution of the minimum critical Reynolds number over the sur- 
face can be easily estimated. The present investigation is limited to 
two-dimensional disturbances. (See ref. 4 for a discussion of three- 
dimensional disturbances . ) 

The distribution of the minimum critical Reynolds number and the 
distribution of the boundary-layer Reynolds number enables the stability 
of the laminar boundary layer with respect to the small-disturbance 
Tollmien-Schlichting type of waves (ref. 4) to be estimated. The 
boundary layer is stable when the boundary-layer Reynolds numbers are 
less than the minimum critical Reynolds numbers and unstable when they 
are greater. If the boundary layer is unstable, the Tollmien-Schlichting 
waves will amplify and eventually cause transition somewhere downstream 
of the location where the boundary layer first becomes unstable. 

It is known that, even though the boundary layer is stable, transi- 
tion can still occur if surface imperfections or other sources of dis- 
turbances generate disturbances sufficiently large to be outside the 
scope of the linear theory (ref. 4) or if the type of disturbances that 
lead to transition are different from those postulated (for example, 
see ref. 5)* Moreover, experiments seem to indicate that extreme cooling 
may cause early transition (ref. 6) although the theory based on the 
Tollmien-Schlichting type of waves predicts that the laminar boundary 
layer on a very cool surface is stable j this phenomenon is not under- 
stood at present. 


SYMBOIS 


A constant 

a velocity of sound 

c wave velocity of disturbance 

_ c 
Ue 


specific heat at constant pressure 



5 


f = ♦ 


m + 1 
2v 0 UeX 


h s “ Cpt + ~ 2 ~ 


ho " OptQ 


m 


Me 


n = 


exponent from Ue 


local Mach number 


= AX m (ref. 2) 

at outer edge of boundary layer. 


^e 


dUe 
dX 

5 0 



correlation number (ref. l) 


R§* minimum critical boundary-layer Reynolds number based on dis- 

placement thickness 6* 

Rq c minimum critical boundary-layer Reynolds number based on 

momentum thickness 9 

h s 

S enthalpy function, — - 1 

ho 

t temperature 



u velocity component parallel to surface 



^e 


— ua 0 

U = , transformed velocity component parallel to surface (ref. 2) 

a e 



k 



U e 


X 

y 



transformed distance along wall (ref. 2) 
distance from wall 


transformed distance from wall (ref. 2) 


Pm 

6 = — y pressure gradient parameter 

m + 1 

7 ratio of specific heats (taken equal to 1.4) 


boundary-layer thickness 


boundary -layer displacement thickness 


, S’-/* 

Jo 


1 - 


pu 


V. w 


n = — l m + — similarity variable 

X 2 v 0 


§ = f *** 1 1 - — jdy, boundary- layer 

JO P e u e \ V 

0^ r = f iL|l - 3-|dY, transformed 

Jo v e \ u e y 


mctrentum thickness 


momentum thickness (ref. l) 


= f f'(l - f' 
J o 


)dT] 


viscosity 


i a 
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V 

p 


Npr 


0 = - 

11 

*fw 

4 


kinematic viscosity 

density 

Prandtl n umb er 



2t'f"M 


♦ stream function (ref. 2) 

Subscripts: 

e at outer edge of boundary layer 

0 stagnation value outside boundary layer 

c at critical layer inside boundary layer, where u = c 


00 


value at which Rq c 


oo when f ' = 1 - — 
Me 


w value at surface 

Primes denote differentiation with respect to tj. Barred quantities 
are dimensional and X, Y are dimensional. 


ANALYSIS 


Derivation of Equations 

In order to calculate the minimum critical Reynolds numbers for the 
similar solutions of references 2 and 3> equations (5.4.3) and (5.4.4) 
of reference 4 are used; these equations can be written as 
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of ref. k) must also be satisfied. It is remarked that the quantity O.76 
in the exponent I.76 in equation (l) follows from the use of a power law 
for the viscosity, with exponent equal to O.76, in the derivation of 
equation (l). 

When the reference length is changed frcm 6 to 0 , equations (l) 
and (2) can be written as 
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and. 


or 



0.58 


( ckA 

-*(*-) c 

^ u bn St 

^y 2 




/ du\ 


\Sy/ 


(4) 


In order to write equations ( 3 ) and (4) in the notation of refer- 
ences 1 and 2 , note that from references 1 and 2 


u 



U 


and 


u c 


a e ^ 
= =£ U e 
®0 


_u_ = _U_ 
u e ^e 


thus 



8 


where 


However, 


(ref. 2) thus, 


Therefore , 


where 


so that 



(5) 


U = U e f ' 

Then 


- i- - i - f (6) 

U e 
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In order “bo obtain the expression for ^u/Sy note that 


8y -,y Sy dy 9 


e 


The definition of 0 in equation (j) is 


r 

1 lsL+ 

1 u e 

If 2 

VqX 


- f - i-W 

J 0 p e Ue \ u e / 


(7) 


or 


= s r°° f(i - f)j- 

p e w o p 0 


ay 


( 8 ) 


but 


± *y = d y 

P 0 a e 

(eq. 6(h) of ref. 2) so that equation (8) becomes 


(9) 


p o a o r 00 

= / f'(l - f ' )dY 

P^ Jn 




but 


fm_+_l Ue 

VqX 


dr) = dY 


2 


( 10 ) 
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Then 


e = 


p_o a o 
p 


e &e _ - r f(i - f)dTi 

+ _1 Ue_ 

2 VqX 

When equations ( 9 ), (10), and (ll) are used, equation (7) becomes 


( 11 ) 


00 

f'(l - f" 

0 


)dri 


( 12 ) 


But 


where 


P _ be 
~ p e 1 


Lli. tijd + s) - 1^1 &' 2 

(which is eq. 31 of ref. 2). Then equation (12) becomes 


dy t 


(13) 




= f f'(l - f 

Jo 


)dq 


where 
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From equation (l4) there is obtained 


a 2 u 




t'f")£a 

by 


where 


With 


it follows that 


by bY by p e 


dg _ A 
by t 

Then 


(15) 



t ' f " ) 


where, from equation ( 13 ) , 



By the use of equation (15) there is obtained 



A 

t 


(16) 


(17) 


(18) 
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When equations (6) and (l4) are used, equation (3) becomes 


R 


.. 1.76 
25Af^t c 


e,c 


bw(^c) 


]f ~ M^(l - 


nr 


f» > 1 _ _ 
c - Mg 

^when M e ^ 1, 


( 19 ) 


where from equation (13) it follows that 


t w - (l + 2L ^ : M|)(l + SJ 


and 


t c ” 


1 + 


4^) 


(1 + s) 


2L ^ : Mi(f , )^ 


( 20 ) 


( 21 ) 


When equations (6), (l4), (l6), and (l'3) are used, equation (4) 
becomes 


-*fw 

■ft 

4 

"3 
f y 


-(tf'" - 2t 


'f") 1 = 0. 


58 


( 22 ) 


The expressions for t, t', ty, and t c in equations (19) and (22) 
Eire given by equations (13)> (l7)> (20), ani (21), respectively. 


Calculation Procedure 

The values of the minimum critical Reynolds number Rq >c were cal- 
culated by means of equations (19) and (22) for the Mach number range 
between 0 and 2.8 for all the solutions witn f^ > 0 presented in 

table I of reference 2 except those for Sy = -1, and for all the solu- 
tions presented in table VI of reference 3* All the solutions of refer- 
ence 3 are for Sy = 0. The special case 3y = -1 is discussed later. 

The values of Rg^ c were also calculated for solutions that are not 
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included in "table I of reference 2 but which are listed in table II of 
reference 2, namely, the solutions for 3 = 0 and = 1, 0, -0.4, 

These solutions were obtained by using the solution for 3=0 
in reference 3 together with the Crocco relation for 3 = 0 , that is, 
s = Sw(l - T')* (See page 3 of ref. 2 .) 

The calculations were made with the aid of the IBM type J 0 b elec- 
tronic data processing machine. Because the value of R e c depends on 

■^c ra i se< i to the fourth power (see eq. ( 19) ) and is thus sensitive to 
the value of f^ and because a high-speed computing machine was avail- 
able, an iterative method was used to find f£. The method was to com- 

the left-hand side of equation (22), for a range of values of q 
beginning with q = 0 . Upon reaching a value of q for which 0 was 
greater than 0.58, this value of 0 and the two preceding values were 
used in a second-order divided-difference interpolation procedure to 
find the value of q at which 0 = O.58. 

Pn a few cases the value 0.58 lay between q = 0 and the first 
value of q in table I of reference 2 ; in these cases two values of 0 
beyond 0.58 were used. Interpolations were made in the tables of given 
data to find f, f 1 , f", S, and S' at this value of q (called q c ). 

The value of f m was also needed (see eq. ( 22 )); this value was obtained 
by the use of equation (l8a) of reference 2 which can be written as 


f»' - pjf '2 - (1 + S)j - ff" ( 23 ) 

The value of Rq^ c was then computed. 

Because near q c , the functions f, f, f", S, and S' are 

usually either monotonically increasing or decreasing whereas the func- 
tion 0 often has a maxi mum and a minimum, the accuracy of the inter- 
polation was improved by using the values of f c , f£, and so forth to 
calculate the value of 0 for q c ; this value of 0 usually differed 
slightly from O.58. A new interpolation to find q c was then made. 

In this interpolation the value of 0 that differed slightly from O.58 
was included in the interpolation and the value of 0 that differed 
most from O.58 was dropped. When the new value of q c was found, inter- 
polations were again made in the tables of given data to find f, f', 
f , S, and S'. A new value of Rq^ c and a new value of 0 were 
then computed. This procedure was continued until 

R 9,c 2 ~ R e,ci 

R e, Cl 


S 0.0001 
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tut never more than six times# Because the data in table X of reference 2 
are given to four significant figures, the final results of the present 
computations were rounded off to four significant figures and are so 
presented in table I. In order to provide "working charts" and to show 
more readily the dependence of Rq c on 3, M^, and Sy, these results 

are also presented in figure 1. 

The case Sy = -1 is a special case because the left-hand side of 
equation (22) cannot be used to compute 0 numerically because the 
quantity ty in the denominator is zero for Sy = -1. (See eq. (20).) 
Equation (22) indicates that, in order that 0 = O .58 when ty = 0, it 
is necessary that either f ’ = 0 or (tf ,M - 2t'f") c = 0. First con- 
sider the condition f^ = 0; the condition (tf m - 2 t’f") c = 0 is dis- 
cussed later. * For M_ < 1, the requirement f £ ^ 1 - -7- does not apply; 

e Me 

thus, any value of f' c between zero and unity is allowable. If the 

quantities that occur in equation (22) are expanded in powers of q and 
only the first power of q is retained, these quantities become 


f' = fyT) 

f" = f^ + fy' T) 

f = -P(l + Sy + S^q) 


t = f 1 + M|j(l + Sy + S^q) 

V - (l * t|)si - 2 

where the result that S() = 0 has been used. (See eq. (l8b) of ref. 2.) 

When these expressions are substituted into equation (22) and powers 
of q greater than the first are neglected, the result for 0, the left- 
hand side of equation (22), is 


0 = 


Jtq 


(1 + Sy ) f 


-fp 

W l- 


3(1 + Sw) + 



(24) 


i-3 CM OJ 'O 



*-1 OJ CVJMD 
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When Sy is sufficiently near -1, the quantities and f^ are 

Q 

both positive and the term 2 S^f^. is much greater than p(l + S^) . 

The quantity 0 therefore increases linearly with q from zero for all 
values of p. As S w approaches -1, the slope of the curves for 0 
against q approaches infinity so that the value 0 = O .58 occurs at 
q = 0. Therefore, q c = 0 and f£ = 0 are allowable values. 

The form of equation ( 19 ) that is valid when q c is near zero is 



If the term p(l + Sv)^ is neglected with respect to 2S^f^ in equa- 
tion (24), a value of O .58 is substituted for 0 and equation (24) is 
solved for q c , the result is 


^c 


0 . 58(1 + Sw) 

2 ^ 


( 26 ) 


If this value of q c is substituted into equation (25), the result is 
an equation for Rg c that is valid for Sy. near -1 and M e < 1} 
namely, 



( 27 ) 


If Sy, is placed equal to -1 in equation (27), the result is that 
Rg^ c = 00 • Thus, for M e 1 and = 0(3^ = -l), the critical 

Reynolds number is infinite. 
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Now consider the condition that (tf ,,, - 2t'f") c = 0. When M e > 1, 

the relation fi > 1 - — must be satisfied. Therefore, q„ cannot be 

^ ,c 

equal to zero and is in fact far from zero for large Mg. The quanti- 
ties f' and t in equation (22) are then not zero. Therefore, in 
order that 0 = O.58 when Mg > 1 and = 0, it is necessary that 

(tf"' - 2t , f") c = 0 ( 28 ) 

The substitution of equation (lj) for t and of equation (17) for t' 
results in a form of equation (28) that contains Mg explicitly, namely, 



Equation (29) can also be written as 



f ltl (f 1 ) 2 - Uf^f") 2 
f"'(l + S) - 2f"S' 


(30) 


When a value of q c is chosen arbitrarily, equation (30) gives the 
value of Mg at which equations (28) and (29) are satisfied. 

Calculations of Mg by means of equation (30) for a range of values 
of q show that equation (28) or (29) is satisfied at two values of q 
for each value of Me above a minimum value that depends on (3 . The 
minimum values of Mg axe found to be greater than unity so that the 
condition (tf m - 2t'f") c = 0 cannot be satisfied for M e < 1. At the 

smaller value of q the relation f ' > 1 - — is not satisfied: at the 

c _ M* 

larger value of q this relation is satisfied when Me is greater than 



^ OJ OJ vo 


3M 
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a value of M e that depends on 3 and is called Mg When Mg is 
greater than M e ^ OT , the larger value of is thus Tj c and is the 

value of t] that is associated with an allowable value of f£, a value 
of f£ for which 


f' > 1 - — 
c = Me 


First consider the case fi > 1 - -i-. Calculations show that for 

Me 

f c > 1 ~ value of M e given by equation (30) increases as f£ 

increases. In order to examine the behavior of Me as f * approaches 1, 
substitute for f ,,f in equation ( 30 ) its expression given by equa- 
tion (23)« Equation (30) then becomes 



As f approaches 1, the quantities f", S, and S 1 all approach zero 
but the quantity f becomes large. Then, considering 


f = 1 

f f = 1 - e 
f" = € 

S = -e 
S' = € 

and keeping only the largest part of each term results in 

^(f , ) 2 [(f , ) 2 - (l + s[| -»3[(f') 2 - l - s] 

fft + 
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3(1 + S)[(f') 2 - (l + sj ->p[(f') 2 - i - s] 
(1 + S)f + 2S' ->f 


Then for f 1 


approaching 1, equation (30) becomes 


+ y ~ 2 ~ “e pjjf ' ) 2 - 1 - £] - ff" 

Lj-± m| 3[(f’) 2 - i - s] - ff" 


(32) 


Thus, as f' approaches unity, the value of Mg that satisfies equa- 
tion (28), or its equivalents equations (29), (30), or (3l)j approaches 
infinity. 

In order to show that the requirement 0 = O.58 is satisfied when 
equation (28) is satisfied and ty = 0, note that the process used to 
obtain equation (32) from equation (3l) shows that the left-hand side 
of equation (28) or (29) is negative for f ' near unity. Then because 
f', t, and f" in equation (22) are positive, the quantity 0 is 
positive for f' near unity. For / 0 the quantity 0 is thus 
zero at a value of f ' and Me given by equation (30) and is positive 
for f' near unity. At the same Mg there is another smaller value 
of f ' at which 0 is also zero but this value of f' is too small 

to satisfy the condition f£ ^ 1 - (See, for example, fig. 2(b) of 

ref. 7.) This smaller value of f* corresponds to the smaller of the 
two values of q mentioned in the discussion that follows the presenta- 
tion of equation (30) • 


By expanding f', f" and so forth around the value of q at which 
0 = 0 and then neglecting terms in q - °f order higher than the 

first, it can be shown by a procedure similar to that used to obtain 
equation ( 2 k) that 0 is approximately proportional to tj - q^^ near 

q = t)0_q. Therefore, because ty appears in the denominator of equa- 
tion (22) the slope of the curve of 0 against f' becomes very large 
as ty becomes very small. Consequently, the value of f ' at which 
0 = O.58 approaches the value of f' at which equation (28) is satis- 
fied. In the limit ty = 0, the quantity 0 is equal to O.58 at this 

value of f . 
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Thus for f^. > i _ A. there is a range of Mg extending to infinity 

for which 0 = 0.58 at t w = 0. At M e = «,, f • = i and tj c = 1. 

Because A, f£, and t c are not zero in this range of Me but ty. is 
zero, equation (19) indicates that R 0 ^ c = ro . Therefore for 

tw = 0(Sw = -1) and a range of Mg that extends to infinity, the 
value of R 0 ^ c is infinite. 

The range of M e determined in this way has a lower limit that 

occurs when f' Q =l--L } this value of Mg is Mg^. In order to find 

Me,®)* equation (30) for Mg must be solved with the condition that 

f c = 1 " I 1 * 16 results of this calculation axe given in table II. 

Note that both conditions that allow equation (22) to be satisfied when 
tw = 0 have been accounted for, namely f£ = 0 when M e < 1 and 

(tf"' " 2t ' f ") c = 0 when M e > 1. (See also page 476 of ref. 7 for a 

discussion of the case tw = 0.) 

For each value of 3 and S*. = -1 there is, in the range of M e 

between unity and the value on the right-hand side of the last column 

of table II, no allowable oscillation in the boundary layer because the 

conditions 0 = O.58 and f' > 1 - -i- cannot be satisfied. The usual 

Me 

interpretation, however, is that the boundary layer is stable in this 
region of M e . Therefore R 0 ^ c =00 for all values of Mg for s* « -1. 

For values of Sw / -l(t w > 0) there can also be a region of M e 
in which there is no allowable oscillation. This region of Mg can be 
found for each value of ^ and 3, when there is such a region, by 
noting that at the upper and lower boundary of the region the conditions 
0 = 0.58 (eq. (22)) and f* = 1 - -L are both satisfied. 

Me 

In order to calculate these boundaries the condition 
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was rewritten as 


1 - M|(l - f*) 2 = 0 (33) 

This term appears in the denominator of equation (19) and, when this 
term is zero, Rq c i6 infinite. Actually, this condition for R 0 ^ c = 00 

is exact and does not depend upon equation (19)- (See page 87 of ref. 4 
and page 469 of ref. 7- ) The calculation was made by choosing a value 
of Mg and then finding f ^ from equation (22). The left-hand side of 

equation (33) was then calculated. This procedure was repeated for a 
range of Mg large enough to allow interpolation for the value of M e 
at which equation (33) is satisfied. This value of Mg is Mg^; values 
of M e>00 are presented in table II and figure 2. The values of Mg j00 
in table II indicate that the upper branch of the curve of Me^oo against 

3 in figure 2 is double-valued between 3 = -0.3884 and 3 = -0.3657 
for Sy = -1 and probably also for part of the range between 3 = -O.3285 
and 3 = -0.3250 for Sy = -0.8. The curve of Mg^ against p in 

figure 2 has been drawn without regard for these double-valued regions. 

It is remarked that, if Mg^ were plotted against fJJ. instead of 3, 

there would be no double values. (See table II of ref. 2 for values of 
f£. ) Note that both conditions that allow R e>c to be equal to infinity 

(eq. 19) have been accounted for; they are 1-y = 0 and 

1 -M 2 (l - fg) 2 =0. The condition fg = 0 occurs together with t w = 0 
for Mg ts 1. 

Figure 3 is a cross plot of figure 2 and shows the connection between 
the wall temperature ratio for R 0 ^ c = 00 when fg = 1 - ~ and M e for 

a range of values of the pressure gradient parameter p. 

Relation Between n, Sy, and 3 

The present results give Rq^ c as a function of the pressure 
gradient parameter 3 and the enthalpy function at the wall Sy. The 
method of reference 1, however, results in a distribution along the body 
surface of the correlation number n which is also a pressure gradient 
parameter but which is not the same as 3. In order to find the dis- 
tribution of Rq c over the surface from the calculated distribution 
of n and the given distribution of S w , it is thus necessary to be 
able to fine 3 when n and Sy are known. 
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In order to find the connection between (3, n, and note that 

(from eq. (22) of ref. l) 


5 2 L 2 _ 


n (l + 3w) 


°tr a U 
^e\^Y 2 , 


Also note that from equation (5) 




ay 2 j v V V Y / * m + 1 


or, upon using equation (10), 


m -f W^2 l.V 


Also note that 


■ I - 

o n TT_ 


=- 1 - e- dY 

0 U e \ U e / 


>. which is eq. (l6) of ref. l) or, upon making use of equation (6), 


p 00 

= / f* 
J o 


(1 - f 1 )cLY 


When equation ( 10 ) is used, this expression for 0 ^ r becomes 


1 n°° 

®tr = " r - / f 1 ( 1 - f * )dT] = 

/ Jo 


m + 1 Ue 
2 - - 


m + 1 u e 
2 z \ 
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Then, equation (34) becomes 


n(l + Sw) = A 2 f"' 


From equation (l8a) of reference 2, it follows that 


( 36 ) 


f w" = "3d + Sy) 


Therefore, equation (36) can be written as 


n = -PA C 


( 37 ) 


The relation (37) was used to calculate r for all the values of 
3 and A given in table II of reference 2; equation (35) shows that 


the quantity A is the same as the quantity 


presented in table II of reference 2. 


^tr\ fm + 1 Ue^ 


which is 


'0 


The relation between n, 3, and Sy is presented in table III and 
figure 


DISCUSSION 


Accuracy 

The values of Rq^ c have been calculated by means of equations (19) 

and (22) which are both approximate. Equation (19) in particular is 
highly approximate and probably is a useful approximation in a range 
of whose upper boundary is only slightly greater than unity. (See 

page 04 of ref. 4.) Moreover, even the more exact method of calculation 
is believed to be adequate only up to a Mach number of about 2. (See 
page 473 of ref. 7*) It is consequently apparent that the present cal- 
culations of R e ^ c cannot be expected to show more than trends with p 

and Sy when M e exceeds unity. 

The accuracy of equation (22) and especially that of equation (19) 
decreases as f,!. increases. The quantity f£ increases as the ratio 
of wall temperature to stagnation temperature increases (Sy increases) 
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and as p decreases, except for cold walls = -0.8). Consequently 
for hot walls and small 3 the present calculations of R 0 c probably 
can only show trends even when Mg is less than unity. 


In order to obtain more direct evidence concerning the 
the calculated values of Re^ c , three comparisons were made. 


accuracy of 
The first 


is shown in figure 5 and is a comparison of the values of R 0 c cal- 
culated by equations ( 19 ) and (22) for the Falkner and Skan profiles 
' ref* 3 ) with the values of R 0 ^ c calculated by Pretsch by an "exact" 

Z th °A (r f- the MaCh number is zero 811(1 th e wall is insulated 
[Me - 0; Sv - 0 ) . The accuracy of the present results is believed to 
be adequate. 


The 


second comparison is shown in figure 6 and is a comparison of 

T. "i n-n rv-P T? ^ t 4- v, > i _ _ *- 


the variation of R 0 , c with Mg for a strong favorable pressure gradient 

and an ( insulated wall (3 = 0.6; S* = 0) calculated by Laurmann (ref. 9) 

by ^oo? XaCt " method with the variation calculated by equations (19) 

+ >/ i* F ? r thiS CaSe the accurac y present calculations seem 

to be adequate up to about M e = 1.3. It is remarked that the theory 

used by Laurmann has been improved by Dunn and Lin. (See ref. 7.) 


The third comparison is the variation with Mg of 


( 'tjA 

\t e / '* 


the ratio 


of wall temperature to the temperature outside the boundary layer required 


for Rq c =00 when f 


c ~ 1 " 7T > vtien the pressure gradient is zero 
AYL e 


( 3 - 0 ). For Mg up to about 2, figure 5-4 of reference 4 shows that 

the variation of with Mg is insensitive to the value of the 

\ e / 00 


Prandtl number and the variation of viscosity with temperature. There- 
fore the accuracy of ( ~ J computed by equation (22) can be tested in 

\ oo 


this range of Mg by comparing these values of (h\ with more accurate 

00 

values even though the Prandtl number is different. Equation (33) is 
also used in the computation of (^\ but is merely a statement of the 
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condition f g = 1 


— . Such a comparison is shown in figure 7j this 

^ . 


figure shows that for M e up to about 2.8 the variation of 



cal- 


culated by the use of equation (22) agrees fairly well with the varia- 
tion given in reference 7* 


It is noted that the indication from figure 7 is that, for values 

/ v 


are too low. The 


of M e greater than about 2.0, the values of I — 

\t e y 

inference is that this result is caused by the use of a Prandtl number 
of unity in the calculations of the velocity and temperature profiles 
of reference 2. This comparison thus shows that formula (22) is adequate 

for the calculation of (^) and M e ,» up to at least Mg = 2. More- 

VeL 

over, the discussions on page 84 of reference 4 and on page 469 of ref- 
erence 7 indicate that formula (22) is much more accurate for the cal- 
culation of ( — | and Mg^oo than is formula (19) for the calculation 

Woo 

of R0 jC * ft is remarked that formula (22) is approximate because the 
number O.58 is used on the right-hand side instead of a function of fg 


and of the velocity and temperature profiles. This function is close 
to O.58 when f£ is small. (See figs. 2 (a), 2 (b), and table 8 of 

ref. 7.) 


The approximate connection between I — - - j and M e is shown in 

VWco 

figure 3 for constant values of the pressure gradient parameter p. An 
increase in p, which means an increase in the favorable pressure gra- 
dient, causes the temperature ratio necessary for Re^ c = 00 to rise 

and also increases the range of Me in which it is possible to make 
r 0 c = «. Figure 3 also indicates that an Lnsulated surface can be 
completely stabilized at Mg equal to about 1.6 if P = 0.4 and for 
a range of M e for P > 0.4. For values of P greater than 0.4, sur- 
faces that are hotter than the insulated surface can also be completely 
stabilized for a range of M e that is centered in the M e region between 

about 1.6 and 2.0 and that decreases as the surface becomes hotter. 
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Because figure 3 is a crossplot of figure 2 , it is not as accurate 
o I? 1 ** 2 ' ^ points of intersection of the curves for 0 constant 

the curves for Sy constant are accurately known but the other por- 
tions of the curves for 0 constant depend on the crossplot. 


Anomalous Results 

The calculations of Rq^ c resulted in two cases in which Rq c 

decreased as 0 increased, an unexpected result. The first case is 
t for Sy - 1 (fig. 1 (a)) when 0 increased from 1.5 to 2 . 0 , a 
large increase in favorable pressure gradient. 

>.4 reas ° n f °r this result seems to be that the length 0 upon 

which Rq^ c is based is sufficiently smaller for 0 = 2.0 than for 

P = 1.5 to cause the decr ease in R 9 , c . Thus, from table II of refer- 
ence 2 , the value of the quantity to which 0 is pro- 

11 v 0 

portional, decreases from 0.1113 at 0 = 1.5 to O.O 6683 at 0 = 2.0 
a decrease of 40 percent. If the reference length had been the dis-’ 
placement thickness, the value of R 5 * >c at Me = 0 would be 14,460 

f ° r v 0 V * 5 W0Uld be l8 ' 290 for P = 2 -°* critical Reynolds 

number Rs*^ c would thus increase with 0 , as expected. 

The second case is that for the highly cooled wall, Sy = -0.8. 

(See fig. 1(d).) In this case R e ^ c decreases with an increase in 0 
for all Mg below Mg^. The two values of 0 that seem to be incon- 
sistent are 0= -0.3285 (f£« 0.0693) and 0= -0.325 (f; = 0.0493). 

This decrease of R e ^ c with increase in 0 has previously been encoun- 
tered and discussed (ref. 10 ) in the comparison between a highly cooled 
two-dimensional, stagnation-point flow ( 0 = 1 ) and a flat-plate flow 

?' with * er ° or smal l rates of mass-flow injection. Note, however, 
that the smallest value of R e , c > that for 0 = 2.0 and Mg = 0 , is 

2.46l x 10^, a value that is larger than any value of R 0 likely to be 

reached. The conclusion therefore seems to be that, for very highly- 
cooled walls with values of R 0 ^ c larger than any value of the boundary- 

layer Reynolds number likely to be met, the effect of a favorable nres- 
sure gradient is destabilizing when Mg < M e> «« Calculations for values 

of Mg up to 8 show that, for values of Mg greater than Mg^, an 

increase in 0 increases R 0 , the usual effect. The values of R fl 

y . c 
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decrease rapidly from R 0 ^ c - 00 to ^0,c ^ ^ or greater than 

M e>00 . (See table I.) 

Because figure l(d), which is for ^ = -0.8, indicates that for 
highly cooled walls at Mg < M e>a8 the critical Reynolds number R e ^ c 

increases as P decreases, the question arises as to what happens as 
the separation point is approached} at the separation point P is neg- 
ative and Rq c is usually near zero. In the solutions of reference 2 
for S„- = -0.&, as the pressure gradient parameter P decreases from 2.0 
to its maximum negative value, - 0 . 3285 , the quantity f w to which the 
skin friction is directly proportional also decreases. A further decrease 
in f", however, is associated with an increase rather than a decrease 
in pT (See table II of ref. 2.) In the region between the value of P 
for separation (f^ = o), namely, -0.3088, and the value - 0 . 3285 , there 
are two positive values of t” for each value of p. Because the skin 
friction is directly proportional to f^-, it is thus f^ rather than P 

which must be used to measure the nearness to separation. Therefore 
R e has been plotted against f£ in figure 8. The two values of P 

that previously seemed to be inconsistent witn the increase in R 0 ^ c 
as P decreases, namely, P = -0.3285 (fw = O.O 693 ) and p = -0.325 
(f|^ = 0.0493) are now seen to be consistent. The conclusion from this 
figure is that, although R 0>c increases as P and f " decrease, a 
value of f" is eventually reached beyond which Re #c decreases rapidly 
with a further decrease in f". The behavior of R 0>c for highly cooled 
walls consequently agrees with the usual behavior, namely, that R 0 , c 
approaches zero as f^ approaches zero at the separation point. 

The data in table I indicate that R 0 ^ c f° r the case P = -0.325, 

S = -0.8 (fJJ. = 0.0^93) behaves in an unusual manner for Me between 
about 1.0133 and 1.116. For Me between 1.0133 and 1.016 the present 
method of confutation results in three values; of R0, c at the same Me- 
(See table I.) The largest values of R 0>c belong to the set that 
increases to infinity at Mg equal to 1.0l6} the other two sets of 
values of R 0 c coalesce at a value of 1,035 X 10 at Me equal to 1.0133* 
If all three values of R 0 ^ c were physically significant, instability 
would occur at the lowest value of R 0jC * Therefore, the physically 
significant value of R 0 ^ c would reach a maximum of 1,248 x 105 at 
M e = 1.0133, decrease dis continuously to 1,035 x 1° this value o± Me, 
and then decrease as shown in table I. Each of the two values of R 0 ^ c 


►q CVJ C\J MD 
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that appears at M e = 1.0133 belongs to a different set of values 
of R 0,c* 0ne set increases with Mg to infinity at Mg equal to 1 . 116 . 

This variation is unlike that encountered for any other case and is prob- 
ably physically unimportant because the values of Ro „ in the other 

Vy C 

set are smaller; this set decreases continuously with Mg in the usual 
manner and is probably the physically significant set. 

CONCLUDING REMARKS 


The minimum critical Reynolds numbers for the similar solutions of 
the compressible laminar boundary layer computed by Cohen and Reshotko 
and also for the Falkner and Skan solutions as recomputed by Smith have 
been calculated by Lin's rapid approximate method for two-dimensional 
disturbances. These results enable the stability of the compressible 
laminar boundary layer with heat transfer and pressure gradient to be 
easily estimated after the behavior of the boundary layer has been com- 
puted by the approximate method of Cohen and Reshotko. 

The previously reported unusual result (NACA Technical Note 4037) 
that a highly cooled stagnation point flow is more unstable than a 
highly cooled flat-plate flow is again encountered. Moreover, this 
result is found to be part of the more general result that a favorable 
pressure gradient is destabilizing for very cool walls when the Mach 
number is less than that for complete stability. The minimum critical 
Reynolds numbers for these wall temperature ratios are, however, all 
larger than any value of the boundary -layer Reynolds number likely to 
be encountered. For Mach numbers greater than those for whi ch complete 
stability occurs a favorable pressure gradient is stabilizing, even for 
very cool walls. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., February 13, 1959* 
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TABLE I.- MINIMUM CRITICAL REYNOLDS NUMBERS FOR SIMILAR SOLUTIONS 
or THE LAMINAR COMPRESSIBLE BOUNDARY LAYER 

(a) S w - 1.0 


1197 

0 

1182 

.2 

11)8 

.4 

1069 

.6 

979-2 

.8 

871.3 

1.0 

749.2 

1.2 

601.7 

1.4 

428.4 

1.6 

258.2 

1.8 

138.4 

2.0 

76.72 

2.2 

47.37 

2.4 

32.84 

2.6 

24.70 

2.8 



P - -0.1 

L 

2.425 

0.7741 

0.1512 

2.431 

• 7758 

.1500 

2.450 

.7812 

.1461 

2.481 

.7097 

.1398 

2.522 

.8007 

.1319 

2-571 

.8134 

.1228 

2.625 

.8267 

.1133 

2.683 

.8406 

.1035 

2.744 

.8543 

.09393 

2.806 

.8673 

.08491 

2.865 

.8790 

•07705 

2.925 

.8901 

.06950 

2.963 

.9002 

.06269 

3.039 

.9092 

.05662 

3.092 

.9172 

.05139 


147-5 

0 

1.692 

0.7273 

O.2727 

142.3 

.2 

1.700 

• 7296 

.2702 

i27.7 

.4 

1.722 

-7371 

.2629 

106 .8 

.6 

1.759 

.7484 

.2516 

84.04 

.8 

1.806 

.7630 

.2370 

63.04 

1.0 

1.863 

•7795 

.2205 

46.20 

1.2 

1.926 

• 7971 

.2029 

34.16 

1.4 

1.992 

.8147 

• 1853 

25.70 

1.6 

2.061 

.8318 

.1682 

19.88 

1.8 

2.129 

.8478 

.1522 

15.84 

2.0 

2.196 

.8624 

.1376 

12.92 

2.2 

2.261 

.8757 

.1243 

IO.76 

2.4 

2.323 

.8876 

.1124 

9.122 

2.6 

2.383 

.8982 

.9075 

.1018 

.09246' 

7*840 

2.8 

2.440 

































TABLE I— MINIMJ4 CRITICAL REYNOLDS NUMBERS FOR SIMILAR SOLUTIONS 
OF THE LAMINAR COMPRESSIBLE BOUNDARY LAYER - Continued 
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TABLE III.- RELATION BETWEEN n, p, A, AND S^. 


2.0 

1.5 

1.0 

• 5 

• 3 

0 

-.10 

-.1295 


Sy — 1.0 

0.06683 
.1113 

.1765 
.2740 
.3334 
.4696 
• 5425 

.5677 


Sy - 0 


n 


-0.008932 

-.OI858 

-.03115 

-.03754 

-.03336 

0 

.02943 

.04174 


2.0 

0.2308 

-O.IO65 

1.6 

.2504 

-.1003 

1.2 

.2761 

-.09148 

1.0 

.2923 

-.08544 

.8 

.3118 

-.07778 

.6 

.3359 

-.06768 

.5 

• 3503 

-.06135 

.4 

.3667 

-.05380 

.3 

.3857 

-.04464 

.2 

.4082 

-.03332 

.1 

.4355 

-.01897 

.05 

.4517 

-.01020 

0 

. 4696 

0 

-.05 

.4905 

.01203 

-.10 

.5150 

.02652 

-.14 

• 5386 

. 04061 

-.16 

.5522 

.04878 

-.18 

.5677 

.05801 

-.19 

.5765 

.06316 

-.195 

.5814 

.06591 

-.1988 

.5854 

.06813 


2.0 

.5 

0 

-.20 

-.24 

-.2483 

-.246 


2.0 

1.5 

.5 

0 

-.14 

-.30 

-.325 

-.3285 

-.3285 

-.325 

-.3088 


Sy = -.4 

0.2944 

.3799 

.4696 

.5544 

.5868 

.6001 

.6045 


Sy = -.8 


0.3551 

.3659 

.4091 

.4696 

.5037 

.5821 

.6107 

.6193 

.6286 

.6335 

.6274 


Sy = -1.0 


-0.1733 

-.07215 

0 

.06i48 

.08263 

.08941 

.08989 


-0.2522 

-.2008 

-.0837 

0 

.03552 

.1016 

.1212 

.1260 

.1298 

.1304 

.1215 


2.0 

0.3833 

-0.2938 

• 5 

.4235 

-.0897 

0 

.4696 

0 

-.14 

.4952 

.03433 

-.30 

.5498 

.09069 

-.36 

.5908 

.1256 

-.3884 

.6400 

.1591 

-.3657 

.6571 

.1579 

-.326 

.6400 

.1335 
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(a) Enthalpy function at the wall, = 1.0. 

Figure 1.- Variation of boundary-layer critical Reynolds number Rq 

c 

with Mach number at edge of boundary layer Me for constant values 
of the pressure gradient parameter 0 . 
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Figure 5.- Dependence of Rq^ c on the pressure gradient parameter 3 

for Mg = 0, = 0. 
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Figure 6.- Variation of critical Reynolds number with Mach number for 
pressure gradient parameter 0 = 0.6 and insulated, surf ace . Sy- = 0. 




b8 
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Figure 8.- Dependence of the minimum critical Reynolds number R 0 on 
the skin-friction parameter for the surface enthalpy parameter 

^ = - 0 . 8 . 
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